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I. INTRODUCTION 



As the classic of multi-scale processes, fully developed Navier-Stokes turbulence remains to be a fortress, which 
besides many heroic attempts has not surrendered to the cohorts of theoretical physicists fiL Hence, any progress 
in a thorough theoretical understanding should be data driven. This means that for a characterisation and better 
understanding of data there is first the need to develop simple heuristic models and to try to learn the most from 
them. Only then, with newly motivated and intriguing ideas and clever theoretical technology, a Trojan horse can 
be found to take this fortress by storm. - This paper is not about the still to be found Trojan horse. It is about the 
revelation of two pitfalls which appear during the early stage of its attempted construction. 

The energy dissipation field in fully developed turbulence reveals striking, close to singular fluctuations. Empirically 
some features of these intermittent fluctuations can be reproduced by multiplicative branching models 0-0]. Their 
particular strength lies in their simplicity and their inherent self-similarity, which proofs useful when discussing 
multifractal behaviour in the limit of a very large Reynolds number 0]. These models rely on two assumptions: 
(al) the weight curdling, associated with the branchings and describing the energy flux from large to small scales, 
is scale- independent, and (a2) different branchings are independent of each other and are not correlated. Of course, 
these assumptions have been directly tested with large Reynolds number atmospheric and wind channel data J7|— [Tc|] 
and conclusions have been drawn that (bl) assumption (al) appears to hold whereas (b2) assumption (a2) appears 
to be in conflict with the experimental findings, where correlations between different, reconstructed branchings have 
been observed. Although some modifications of the cascade models have been proposed to include these effects [0, 
many physicists think of it as a major drawback to the cascade models. 

This criticism of the cascade models is by no means justified and this brings us back to the pitfall-business mentioned 
in the beginning. The immediate testing of the cascade model assumptions (al) and (a2) from data is not as straight- 
forward as had been long anticipated. In order to explain this point we have to understand how the measurements 
are done: 

An anemometer records one component, which is usually longitudinal, of the velocity field as a time series in one 
point. Given that the fluctuation strength is small this time series can be interpreted as a one-dimensional spatial 
realisation of the velocity field at a fixed time; this is known as Taylors frozen flow hypothesis. Since the energy 
dissipation tensor is a functional of spatial derivatives of the velocity field, the spatial realisation of the streamwise 
component of the velocity field can be converted into a spatial realisation of the streamwise components of the 
energy dissipation tensor. It is the statistical properties of this one-dimensional energy dissipation field which are 
compared to, for example, the multiplicative cascade models. The important thing to note here is that the energy 
dissipation field is measured one-dimensional, but is generated from the three-dimensional dynamics of fully developed 
turbulence. If we assume a cascade dynamics to hold in three dimensions, then the energy flux from large to small 
scales is conserved in every three-dimensional branching. For a one-dimensional cut through the three-dimensional 
dynamics this, of course, needs not to hold. As a consequence, branchings used in one-dimensional multiplicative 
cascade models should not conserve the energy flux locally. Globally, that means on average, conservation of energy 
flux should be fulfilled. This then means, that the weights (or dynamical multipliers) used in the cascade evolution 
can not be reconstructed from the one-dimensional cut through the energy dissipation field resolved at the finest scale. 
With other words, the reconstructed multipliers obtained by smoothing from fine to large scales are not equal to the 
dynamical multipliers of the cascade evolution from large to fine scales. 

A second pitfall occurring in the comparison of cascade models with data has already been noticed previously 
fj"lf : multiplicative cascade models lead to hierarchical structures which are not homogeneous. The time series 
measurements sketched above are unable to trigger on these hierarchical structures. Hence, homogeneity has to be 
restored in the cascade models first, before they are compared to data. 

For a fair comparison between cascade models and data these two points have to be taken into account: (cl) 
no local conservation of energy flux in the branchings of one-dimensional multiplicative cascade models since they 
represent one-dimensional sections through an underlying three-dimensional cascade dynamics, and (c2) restoration 
of homogeneity in the hierarchical multiplicative cascade models. Taking these two points into consideration we will 
show in this Paper that multiplicative cascade models can be constructed, which still respect assumptions (al) and 
(a2) and which at the same time are successful in explaining the presently available experimental observations (bl) 
and (b2). 

The organisation of the Paper is as follows: In Sect. II we summarise the experimental facts about multiplicative 
cascade processes in fully developed turbulence. Section III already reveals that three-dimensional multiplicative 
cascade models, observed in one dimension, lead to fix-point solutions for the various reconstructed multiplier distri- 
butions which come close to the experimentally observed ones. Some analytical insight into the fix-point behaviour is 
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given in Sect. IV. Sect. V considers three specific popular variants of one-dimensional multiplicative cascade models, 
the asymmetric binomial, the log-normal and the log-Poisson model, and studies the implications of our two points 
(cl) and (c2). Conclusions are given in Sect. VI. 



II. EXPERIMENTAL FACTS ABOUT TURBULENT CASCADE PROCESSES 



A. Discrete multiplicative branching processes 



Multiplicative branching processes serve to model properties of the intermittent fluctuations observed in the energy 
dissipation field of fully developed turbulence. They relate the energy flux El at some integral scale L to E r = 
ElWi(X) ■ ■ ■ Wj(X) contained in a subinterval of size r — X~ J at scale J by a product of mutually independent 
random weights Wj (A). The length ratio A represents a free parameter; the most simple and most commonly chosen 
value is given by A = 2. Adopting also this latter choice and dropping the dependence of the random weights Wj on 
A, the multiplicative branching processes are from now on denoted as binary multiplicative cascade processes. 

More precisely now, the energy flux Er' , contained in the interval < k < 2 3 with length r = L/2- 7 splits into a left 
(L) and right (R) offspring interval, each of length r/2, whereby the content propagates according to E^ 1 ^ = WlE^ 
and -Ejfe+i^ = WrE^ , respectively. For each breakup the two random weights Wl,r > are chosen, independently 
from any preceding breakup, according to a joint probability density p(Wl, Wr). The splitting function, as we denote 
the latter, guarantees a full description of a binary breakup. For the special case where the splitting function is 
concentrated along the diagonal Wl + Wr = 1, i.e. p(Wl,Wr) = p(Wl)S(Wl + Wr — 1), each breakup strictly 
conserves energy flux. For all other forms the relation E^ 1 ^ + ^fc+i = nee ds not to hold; however, on average 
it should hold. 

A particular simple example is the binomial splitting function 



P(W L ,W R ) = ~ 



S(Wl-^)+S(W l -^—^) 



5{W L + Wr ~ 1) , (1) 



which is known as the p-model Only two different local breakups are possible, each with the same probability. 
Here, energy flux is conserved for every local splitting. 

In general, multiplicative cascade models are based on two assumptions: (al) the existence of a scale-independent 
splitting function p(Wl,Wr) and (a2) statistical independence of the random weights Wl,r at one breakup from 
those of any other breakup. Once p(Wl, Wr) is chosen, the assumptions (al) and (a2) allow to determine all moments 
and scaling exponents of the energy dissipation e r —E r /r by inverse Laplace transforms |l^ , [l3|] . Many experimental 
measurements of moments and scaling exponents confirm that this simple construction reproduces the measured 
multifractal aspects of the energy dissipation field amazingly well (e.g. ||). 



B. Experimental facts 

From a theoretical point of view the various parametrisations suggested for the splitting function p(Wl, Wr) serve 
to reproduce the scaling exponents observed in the one-dimensional cuts of the energy dissipation field, the latter being 
converted from velocity time series obtained from anemometers and employment of Taylors frozen flow hypothesis. 
Within the experimental error bars, such different models as the p- model ||, the log-normal model jl4|,[l5), the log- 
Poisson model [p~6| |l~S|] and others are all capable to describe the observed scaling exponents. As a consequence of 
this model insensitivity some efforts have been taken to extract information about the splitting function directly from 
data 0-^]. In the following we will briefly sketch the operational procedure and summarise the essential results. 

In the inertial range (0 < j < J) the energy dissipation e(x, r) at length scale r = LjV and position x = kL/2^ is 
defined by 

= * (* = = |) = I ^ s(x')d X ' , (2) 
where < k < 2 3 . The experimentally measured 'backward' energies then become 
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(fc+l)2 J -J'-l 

4 j) = §^= E E \ J) • (3) 

They are obtained by successive summations from the finest resolution scale J to larger ones and are generally not 
equal to the 'forward' energies E£ J , which arise as intermediate states in the evolution of the cascade from larger to 
smaller scales. For a clearer distinction we denote the former with a bar. - Based on (||), so-called left, right and 
centred multipliers at scale j and position k are now operationally defined as 

They may be denoted as base 2 multipliers, because the associated change of interval lengths refers to A = 2. They 
depend on the relative position of parent to offspring intervals. Multiplier distributions p{M^>) are obtained by 
histogramming all < k < 2 J possible multipliers at a given scale and averaging over many realisations. 

For cascade models with an energy flux conserving splitting function p(Wl,Wr) = p(Wl)S(Wl + Wr — 1) the 
random weights Wl/r are directly related to the multipliers M^/r. Due to local energy flux conservation we have 

E 2k 1] + E 2k+i = ( w kl + w i%) E k ] = E k ] = E k \ so that the 'forward' energies E [ k j) are then identical to the 
'backward' energies E { k j) . Hence, = {W^ L E^)/E^ = wjfl, and, analogously, M kR = W^ R . For this case the 

multiplier distributions p(M^J R ) reproduce the scale-invariant splitting function p(Wl, Wr). 

Experimental analyses performed in atmospheric ]7],|| and wind tunnel turbulence Q have revealed that within the 
inertial regime the multiplier distributions p(M^j R ) are to a very good degree independent of the length scale. The 
following parametrisation has been found empirically |g] : 

P(M L/R ) = 1 ^[M L/R (l-M L/R )f- 1 , (5) 

where T stands for the gamma function and where (3 « 3.2. See Fig. 1. In order to link this scale-invariant multiplier 
distribution directly to a splitting function, local energy flux conservation in each binary breakup of the cascade 
process has to be assumed. This then leads to 

p{W L , Wr) = [WlWr] 13 ' 1 5 (W L +W R -1) (6) 

and the assumption (al) of Sect. II. A of the existence of a scale-independent splitting function for multiplicative cas- 
cade models appears to be confirmed. Notice, however, that the additional assumption of local energy flux conservation 
had to be employed to convert the multiplier distribution (^|) into the splitting function (j^). 

What about assumption (a2) of Sect. II. A? If the multipliers M k \/ R are independent from those at other scales 
and positions, as required from the multiplicative cascade models, then the conditional multiplier distributions 
P(^2kt+i) L/R^k \/b)i wmc h correlate an offspring multiplier with its parent multiplier, should not depend on the 
latter, i.e. K^2fc(+i) L/R\^kL/n) ~ K^ifcf+i) l/r)- ^ ne experimental findings, however, indicate a dependency (see 
Fig. 1): given a small parent multiplier M^ L , R < 1/2 the conditional multiplier distribution p{M^^ \M^J R ) is more 
narrow than the unconditional multiplier distribution p{M^ R Vl ), and, given a large parent multiplier M^' L , R > 1/2 

p(M^J R V> \M^J R ) becomes broader than p(M^ R ^); here averages have been taken over all four possible combinations 
LL, LR, RL, RR of the parent and offspring multipliers. - This observation has led to the conclusion about turbulent 
cascade processes, that assumption (al) appears to hold, whereas assumption (a2) appears to be in conflict with 
the experimental findings. To save some features of multiplicative cascade models further extensions such as the 
correlated p-model j^] have been suggested, which attempt to account for the effects seen in the conditional multiplier 
distributions by introducing a dependency of the local cascade branching on the magnitude of the local strain rate. 

From an experimental point of view a multiplier analysis needs not only be focused on left- (or right-) sided 
offsprings of a parent interval; any relative position is qualified. Such a generalisation has been proposed by Novikov 
|fl^| . In this context centred multipliers are most often considered ^,0; in our notation their definition is given in 
Eq. (|J). It has been observed that the unconditional centred multiplier distribution p(M^) is also scale-independent 
in the inertial regime, but is more narrow than the unconditional L/R- multiplier distribution p(M^J R ); see Fig. 
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2. The latter effect reflects the non-homogeneity of the breakup, ft can be explained, at least in a qualitative 



(j+2) 



manner, within the binary multiplicative cascade processes: writing E^ k+l 



W 2k,R h 2k 



AS) 
'k 



and E 



(J+2) 



J 4fc+2 ~~ ^ 2fc+l,L JJ 2fc+l — "2fe-|-l,£ 



— WoiT^r-Eol^V = Wo^ti r wjf R E k J> , respectively, and employing the energy flux-conserving splitting 



it?) 
'k 



), the distribution p(M k J ^,) of Eq. (^) becomes 



p(M$) = 



W U+1) ) 

vy 2k,R I 



k,C 



(i+i) WU+ 1 ) \ 

2fc+l,£' Vv 2k+1,R> 



w. 



0+1) wGO 



2/. J,' 



fe,L 



1,L vv k.R 



2k- 



(7) 



The result is shown in Fig. 3. 

In addition to the unconditional centred multiplier distributions also a conditioning on centred parent multipliers 
has been investigated || . A small centred parent multiplier leads to a narrowing of the unconditional distribution; 
in addition the distribution is not anymore symmetric and is skewed towards smaller multipliers. A large parent 
multiplier, on the other hand, leads to a broadening of the distribution with an asymmetry towards larger multipliers. 
Consult again Fig. 2. 



III. THREE DIMENSIONAL BINOMIAL CASCADE MODEL AS OBSERVED IN A 

ONE-DIMENSIONAL WORLD 

The comparison between the multiplicative character of binary branching processes and the experimental multiplier 
analysis is not as straightforward and clear as it has been summarised in Sect. II. B. In fact, this comparison is not 
justified in a rigorous way as it neglects two important features: (cl) violation of energy flux conservation in the 
bivariate splitting function and (c2) translational invariance. While the latter is the subject of the second Subsection, 
we now concentrate on the former. 



A. Consequences of local non-conservation of energy flux 

In order to identify multipliers Ml /r with weights Wl /r of the bivariate splitting function, energy flux conservation 
(Wl + Wr = 1) of the latter had to be assumed. This assumption is not well justified. The experimental analysis is 
restricted to measure the energy dissipation field of three-dimensional turbulence only along a one-dimensional cut. 
For the modelling this implies that, for example, if we were to consider a three-dimensional multiplicative branching 
process with energy flux conservation in every local three-dimensional breakup, then in one dimension the same 
process would appear as to violate local conservation of energy flux. 

Take, for example, a three-dimensional binomial model, where a cube splits into eight subcubes and where four 
randomly chosen subcubes get a large weight W+ = (1 + a)/8 and the other four get a small weight W- = (1 — a)/8. 
Along a one-dimensional cut this three-dimensional branching is described by the bivariate splitting function 



p(w l ,w r ) = Us(w l -^ 



a 



5 W L - 



1 — a 



p{W L ) 



^p(Wr) + ^S{W l + Wr-1) 



4 / 1 - a 



1 — a 



(8) 



where 



P{W L/R ) = ^ 



W L/R - 



+ S W L/R - 



1 — a 



(9) 



The second contribution in (g) represents a small anticorrelation between weights Wl and Wr and is reminiscent 
of energy-flux conservation in three dimensions. However, the uncorrelated, first contribution in (|^) dominates and 
clearly does not respect conservation of energy flux in one dimension. 
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Note, that for a one-dimensional binary multiplicative cascade process with the splitting function (g) the forward 
and backward energies E ( / } and E<jp are no longer identical. This is illustrated in Fig. 4, where the (forward) 
evolution of one cascade realization down to scale J is compared with the (backward) averaged field. Only at the 
finest resolution scale the forward and backward fields are identical. It is due to the non-conservation of energy flux 
in the splitting function that the two fields are different for all other resolutions < j < J. Now, the forward field 
is faithfully described by the splitting function p(Wl, Wr) depending on the left and right weight Wl and Wr. For 
the backward field it is the left and right multipliers Ml and Mr which are being deduced. Since the two fields are 
now different, the multiplier distributions p(M l /r) have no longer a simple relationship with the splitting function 
describing the physical branching of the cascade; the former will be different from the latter. This can be seen in 
Fig. 5a, where an averaging over 10 5 simulated cascade realisations has been performed. While the weights Wl/r 
in the forward evolution with the splitting function (||) may take only two values (1 ± a)/2, the multipliers M L / R 
defined in Eq. (^]) take more and more distinct values as the scale difference J — j increases and becomes soon quasi- 
continuous. Already after three backward steps the left (and right) multiplier distributions apparently converge to a 
scale-independent limiting form which comes surprisingly close to a parametrisation of the experimentally observed 
multiplier distribution shown as continuous line. 

Since the splitting function (^) is symmetric under exchange of Wl «-> Wr, the unconditional right multiplier 
distributions are identical to the left ones illustrated in Fig. 5a. Of course, this does not hold anymore for the centred 
multipliers Mc', their unconditional distributions are depicted in Fig. 5b. As J — j increases these distributions also 
converge to a stable limiting form, which is now more narrow than the limiting form experienced for the unconditional 
left (or right) multiplier distributions. At least, and this is very encouraging for the moment, the trend is right when 
we think of the experimental facts presented in Sect. II. b. 

The scale-independence of the limiting forms of the various unconditional multiplier distributions suggest that they 
can be understood as solutions of a fix-point problem. In Sect. IV we will further elaborate on this suggestion. For the 
moment we are further driven by curiosity and ask whether we already observe some effects in conditional multiplier 
distributions. The simulations reveal that those themselves converge to limiting forms after only three or four backward 
steps, starting from the finest resolution scale J. Figures 6a and 6b show conditional multiplier distributions at scale 
j = 3 < J = 9. The multiplier distribution p(Mj?J R \Mj?ZP) conditioned on the previous multiplier M L ^ R \ where an 
average over all four possible daughter/parent combinations is performed, reveals a dependence on the latter. For a 
small parent multiplier, here 0.2 < Mj?,^ < 0.4, the distribution becomes somewhat broader than the unconditional 

(j—X) 

one, and for a large parent multiplier, here 0.6 < M L , R < 0.8, the conditional distribution becomes a little more 
narrow. Although, in view of the experimental findings, these effects are in the right direction, we resist not to 
overstress them, since the associated magnitude in the deviations from the unconditional distributions appears still 
too small. A similar result is obtained for the centred multiplier distribution p{M < ^p\M < ^~ 1 ^ 1 ) conditioned on the 
parent multiplier M^~ x \ For Mg 1 ^ small, here < Mq 1 ' < 0.5, the distribution broadens, and for a large parent 
multiplier, here 0.5 < Mq < 1, it narrows to some extend; note, however, that the experimentally observed shifts 
to the right and left, respectively, are missing. 



B. Consequences of translational invariance 

Due to their hierarchical structure discrete cascade processes are obviously non-homogeneous, n-point correlation 
functions of the generated energy dissipation field are not translationally invariant; see Ref. for a visualisation. 
This is in conflict with the experimental findings, where n-point correlation functions are deduced from velocity time 
series obtained from anemometers and employment of Taylors frozen flow hypothesis. Experimental realisations of the 
energy dissipation field are obtained by chopping, more or less arbitrarily, strings with length approximately equal to 
the correlation length L out of the recorded time series. For cascade model realisations, however, we have knowledge 
about where the realisation of length L begins and ends. Hence, in order not to compare "apples with oranges" , 
we have to introduce random translations within model realisations. Evidently, this will also influence multiplier 
distributions. 

Two slightly different schemes have been used so far to restore homogeneity in the (discrete) cascade models. In 
Ref. fll[] two independent cascade realisations, each having been evolved J cascade steps, have been attached and 
an observation window, having the same length as one cascade realisation, has been randomly placed within the two 
attached realisations. Sampling over many realisations and many random shifts then leads to homogenous n-point 
correlation functions. Moreover, it has been found that the multiplier distribution (|s|) associated with the splitting 
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function (|6j) remains invariant under this scheme. The same findings also hold for a slightly modified scheme, which 
has been suggested in Ref. jL9| and which we will employ in the following: for the target resolution scale J of a given 
integral length scale L a longer cascade realization with J + 3 steps is generated (corresponding to an integral length 
scale 8L) and an observation window of size L is shifted randomly by t bins within 8L. In other words, only the 
2 J bins of the generated E^ +3 \ which lie within the randomly placed observation window are considered, giving 

a simulated and translated — sifj^ 3 , where t is a uniformly distributed integer within [0, 7-2 J — 1]. Then the 
multipliers are determined again by (|[) and (|J) and sampled over a large number of configurations and random shifts 
t. In the simulations the second scheme to restore homogeneity will be applied with J + 3 = 12 cascade steps, and 
the distributions have been averaged over 10 5 independent realizations. 

The results for the various multiplier distributions, determined by cascade simulations with the splitting function 
(|J) and with the adopted scheme to restore translational invariance, are presented in Figs. 7 and 8. For intermediate 
scales, here 1 < j < 5, the unconditional left multiplier distribution p(M^) converges again to a stable limiting 
form, which is now in nearly perfect agreement with the experimentally deduced left multiplier distribution (||) with 
(3 = 3.2. Apparently, the applied scheme to restore homogeneity acts as to narrow the distribution to some extend; 
compare Figs. 5a and 7a. Note a small detail: for very small scales, here j = 0, the actual multiplier distribution starts 
to deviate again from the limiting form experienced for intermediate scales; this narrowing effect has already been 
pointed out in Ref. pT] , where a slightly different scheme to restore translational invariance had been applied to a 
binary cascade with the splitting function (^). The differences between conditional left/right multiplier distributions, 
which are conditioned on different values of parent multipliers, arc about the same with and without the proposed 
scheme to restore translational invariance; confer Figs. 6a and 8a. The unconditional centred multiplier distributions 
also show convergence to a slightly modified limiting form, but instead of a further narrowing they experience a 
small, but noticeable asymmetry towards smaller multipliers; compare Figs. 5b and 7b. This is a consequence of the 
associated conditional centred multiplier distributions (Fig. 8b): for a small centred parent multiplier the conditional 
distribution narrows as in Fig. 6b, but this time it is also skewed towards smaller centred daughter multipliers; a large 
centred parent multiplier leads to a broadening of the conditional multiplier distribution with a simultaneous shift 
towards larger centred daughter multipliers. 

Overall, the various multiplier distributions become more realistic once restoration of translational invariance is 
accounted for in the specific discrete cascade model with splitting function (||), which is motivated from a three- 
dimensional binomial cascade process as observed in one dimension. Even the effects seen in the conditional multi- 
plier distributions go, at least, in the right direction. Evidently, a non-energy flux conserving splitting function and 
restoration of homogeneity are two key inputs to explain the current experimentally observed multiplier phenomenol- 
ogy within the simple cascade models. A further key input, namely asymmetric splitting functions, will be discussed 
in Sect. V. But before that, we intend to enlighten the revealed fix-point behaviour in the multiplier distributions and 
discuss a cascade, where the splitting function is given by a symmetric Beta-distribution. 



IV. TURBULENT CASCADE MODEL (STUDY I): BETA WEIGHT DISTRIBUTION 

The one-dimensional splitting function (^|) of a three-dimensional binomial cascade process contains an anticorre- 
lation between the left and right weight. It is the result of local conservation of energy flux in three dimensions. 
However, this anticorrelation remains small and, when compared to the leading term in the one-dimensional splitting 
function, may well be neglected. The one-dimensional bivariate splitting function then takes on a factorised form: 

p(W l ,Wr) = P (Wl)p(Wr) . (10) 

From now on, this means Sects. IV and V, we will only discuss splitting functions belonging to this class. In this 
Section we begin with a Beta-distribution for p(Wl/ r ): 

P(W L/R ) = ^^[W L/R {1-W L/R )] - 1 . (11) 

One motivation for this distribution comes from the experimentally observed left (or right) multiplier distribution, 
for which a Beta-distribution with (3 = 3.2 represents a good parametrisation; see discussion around Eq. (|^). Another 
motivation will be outlined in Subsect. IV. b in connection with a first analytical insight into the fix-point behaviour 
of multiplier distributions. 



7 



A. Multiplier distributions without translational invariance 

Again, from simulations we extract the various multiplier distributions belonging to a binary discrete cascade model 
with the factorised splitting function (|ll]). In this Subsection we do not employ a scheme to restore translational 
invariance; this will be done in Subsect. IV. C. Both, the unconditional left as well as centred multiplier distributions, 
illustrated in Fig. 9, converge to respective limiting forms, where the latter distribution is more narrow than the 
former. Only at the very finest scales, which means a large scale index j, the distributions deviate from the limiting 
forms. At first thought, it might look puzzling why the Beta-distribution (/3 = 3.2), which is used for the splitting 
function, is not reproduced for all scales in the unconditional L/R- multiplier distributions. Note, however, that the 
used splitting function has a factorised form, so that there is an enhanced probability that the two weights Wl and 
Wr are nearly equal; since then Ml w Mr ps 0.5 holds for the multipliers, this explains the narrow L/R- multiplier 
distributions for j close to J, the total number of performed cascade steps. If we had chosen a splitting function 
like @, reflecting a Beta-distribution with local conservation of energy flux, then the extracted unconditional L/R- 
multiplier distributions would have been the same Beta-distribution for all scales j. On second thought we should 
then ask, why does the L/R- multiplier distribution converge back to the Beta-distribution? In the next Subsection 
we will work towards an answer of this question. 

The conditional multiplier distributions p(M^j R \M^ R 1 ' > ) and p{Mq\Mq ^) also show convergence to limiting 
forms. However, these limiting forms seem to be independent of the respective parent multiplier; this is illustrated 
in Fig. 10. This outcome differs from the results obtained in the previous Section. Whether effects in conditional 
multiplier distributions occur or not, apparently seems to depend on some specific features of the chosen splitting 
function. 



B. Fix-point problem 

In the previous Subsection as well as previous Section we have started to realize that various multiplier distributions, 
resulting from one-dimensional cascade processes derespecting local conservation of energy flux, appear to converge 
to scale- independent limiting forms. These limiting multiplier distributions may well be understood as solutions of 
an associated fix-point problem. We will now formulate this fix-point problem and discuss it for the particular case 
of a factorised Beta-distribution splitting function. 

We focus on the unconditional left multiplier distribution. We also do not consider the scheme to restore translational 
invariance, since it renders an analytical approach to be impracticable. According to Eq. (||) the left multiplier L 
at scale j is defined as the ratio of two resumed energies E^ 1 ^ and Er^ . The resumed energy Er^ is equal to the 

product of j weights, drawn during the first j cascade steps, with the total resumed energy E^ J ~^ of a subcascade 
with only J — j cascade steps: 



E^=W?) R ...W^El J -» . (12) 



For the other resumed energy we find 



Eg +1) =wf] R ...W^W^Eg-^ . (13) 



The first j weights W^j R , W^, R in (12) and (|T^ ) are identical. The total resumed energy E^ J can be 



expressed as a weighted sum of E^ k J ' X ' of the left subbranch and Ej^A l ' of the right subbranch: 

E^=W^E^ + W^E^ ; (14) 
this follows from (||). Note, that E^ - - 7-1 -* = E^T^ and E R ~^^ = E^+i are independent variables. Com- 



bining ( |l^ ) - (0), the left multiplier then becomes: 



M 



' ~ W^E^-^+W^E^-V 



(15) 



8 



The quantities wjf\ W R , E^~*~^ , ^~ X ' are independent statistical variables; the first two of them are drawn 
from a scale- independent splitting function p(Wl, Wr) = p(Wl)p(Wr) . Since within the simulations the limiting form 
of the left multiplier distribution appears to be scale-independent, this suggests that the statistics of the resumed 
energies El/r should also be scale-independent. 



Denoting the probability density function of the resumed energy as p{E), we deduce from (14) 

p(E (f) ) = J dW L dW RP (W L )p(W R ) J dE { /- 1) dE^'- 1) p(E { /- 1) )p(E^'- iy ) 

6 (eU"> - [W L E ( /~ X) + WnEg-V] ) , (16) 

where the substitution j' = J — j has been introduced. A Laplace transformation translates this equation into an 
equation for the corresponding characteristic function: 

Z^[X] = (exp(\EV'ty = J p(E^">)e x&J '\lE^'^ = (^J dWp(W)Z^'-^[WX}^ . (17) 

If, here, we could find a scale-independent fix-point solution, then this would explain the scale-independent limiting 
form of the left multiplier distribution. 

For the case of a splitting function with local conservation of energy flux, i.e. p(Wl, Wr) = p{Wl)S(Wl + Wr — 1), 
it is simple to find the fix-point solution. Eq. ([Tt|) then generalises to 



Z^[X] = I dW L dWR P {W L ,WR)zV'-V[W L \} zW-V[WrX] 

= [ dWp{W)Z^'- 1 \W\]Z^'-^[{l- W)\] , (18) 



which is fulfilled by ZU">[\] = Z[X] — expA. Of course, then p(E) = 5(E — 1) and the left multiplier (|l5| ) becomes 
Ml = Wl, so that its distribution is scale- independent. We know this simple result already from before, because for 
an energy flux conserving splitting function the left/right multiplier is always equal to the left/right weight. 

For the case of a factorised splitting function, which of course does not respect local conservation of energy flux, 
the fix-point solution of (|l7]) is much more difficult to find. In fact, for the moment we are not aware of a general 
solution. We only know about a solution for a specific choice of the splitting function [^0| : the Beta-distribution given 
in Eq. ([ll]). The fix-point characteristic function Z[X] for the resumed energy is then the characteristic function of a 
Gamma-distribution: 

)>-'» 



As has been demonstrated in Ref. |2Q| , this can be checked by inserting ([ll]) and ( |l9|) into(0). Just to be on the safe 
side, we have also performed a cascade simulation with j' 3> 1 steps to determine the probability density function of 
E^ ); the result is in perfect agreement with the Gamma distribution 

P(E) = ^E^e Xp (-2(3E) (20) 

associated to the characteristic functio n (|l9| ) . - According to ( |l5| ) the convolution of the Beta-distribution ([ll]) for Wl 
and Wr with the Gamma-distribution (GO) for El and Er leads to the limiting form of the left multiplier distribution. 
As expected this convolution matches the limiting form (of Fig. 9) obtained from cascade simulations. 

The fix-point distribution p{E) for the resumed energy not only explains the scale-independent unconditional left 
(or right) multiplier distribution, it also leads to scale-independent conditional multiplier distributions. Changing 
notation as explained in Fig. 11, a left daughter and parent multiplier become: 

Mll = Mf = WllEll 

WllEll + W lr E L r 

Ml = M«-V = WlWllEll + WlWlrElr 

L WlWllEll + WlWlrElr + WrWrlErl + WrWrrErr 
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All six weights W L , Wr, W ll , Wlr, W rl , W rr as well as all four resumed energies E LL , E LR , Erl, E rr are 
independent variables and are drawn from scale-independent probability density functions p(W) andp(_E), respectively. 
Hence, the probability density p(Mll, Ml) is also scale-independent. Since the two multipliers Mll and Ml partly 
consist of identical random variables, they are correlated. This means that p(Mll, Ml) ^ p(Mll)p(Ml), which is 
responsible for the correlation effects seen in the conditional left multiplier distribution p(Mll\Ml); however, as has 
been realized in Fig. 10a, these effects remain unnoticeably small for the special case of the Beta-distribution ( p"l| ) 
entering into the splitting function (|l0|). - Note again, for an energy flux conserving splitting function we would 
always have Ml = Wl and Mll — Wll, so that p(Mll, Ml) — p(Mll)p(Ml) with no correlations. 

With similar arguments the limiting forms of the unconditional as well as conditional centred multiplier distributions 
can be traced back to the fix-point probability distribution p(E) of the resumed energy and the scale-independent 
splitting function p(Wl, W r ). For example, in view of Fig. 11 the centred multiplier is expressed as 

M = WlWlrElh + W r W R lE R l 

C WlWllEll + WlWlrElr + WrWrlErl + WrWrrErr ' 

which clearly reflects a scale-independent distribution. 



C. Multiplier distributions with translational invariance 

Now, that we have gained some analytical insight into the occurrence of scale-independent limiting forms for the 
multiplier distributions, we again rely on simulations and study the effect of restoring translational invariance onto 
various multiplier distributions resulting from the factorised splitting function ( |l0|) with (|ll|). As in Fig. 9, Fig. 
12 illustrates the unconditional left and centred multiplier distribution. The scale-independent limiting form of the 
unconditional left multiplier distribution becomes again slightly more narrow after application of the scheme to restore 
translational invariance, whereas the corresponding centred multiplier distribution remains practically unchanged. 

The limiting forms of the conditional left /right multiplier distributions show practically no dependence on the 
parent multiplier; only for extreme small values of the latter some correlation effects are observed. See Fig. 13 for 
an illustration. For the conditional centred multiplier distributions the scheme to restore translational invariance has 
some more impact: the daughter and parent multiplier become positively correlated, which means that for large/small 
parent multipliers the occurrence of large/small daughter multipliers is enhanced. However, no broadening/narrowing 
of these distributions is observed. 

We conclude that a Beta-function parametrisation of the factorised splitting function is not able to describe all 
effects observed in various multiplier distributions. However, for this special case we succeeded to understand the 
fix-point behaviour. In the next Section we will study some more realistic splitting functions. 



V. TURBULENT CASCADE MODEL (STUDY II): ASYMMETRIC WEIGHT DISTRIBUTIONS 

Adopting the factorised form p(Wl, Wr) = p(Wl)p(Wr) for the splitting function there is more to expect for the 
weight distribution p(Wl/r) from the one-dimensional section through a three-dimensional energy flux conserving 
cascade process. In an extreme three-dimensional branching event the parent energy is completely transfered into one 
subcube; the other seven subcubes getting nothing. The energy density contained in this one subcube is now eight 
times as large as before, so that the energy contained in the projected one-dimensional subinterval is four times as 
large as in the one-dimensional parent interval. Hence, the values of the weights Wl/r should not be restricted to 
< Wl/r < 1, but to < Wl/r. < 4. Since the average weight will remain (Wl/r) = 1/2, the weight distribution 
p(Wl/r) should then become asymmetric around 1/2. - In this Section three factorised splitting functions, having 
different asymmetries, are picked and their implications on the various multiplier distributions are discussed. We only 
show results obtained after application of the scheme to restore translational invariance. 

A. Asymmetric binomial weight distribution 

The simplest example of an asymmetric splitting function p(Wl, Wr) = p(Wl)p(Wr) is a modification of (jl|), (||) 
and (||): 
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P(W L/R ) = -%-8 (w L/R - i-^) + -J5_* (w L/R - !±f») , (23) 
ai + a 2 V 2 / ai + a 2 V ^ / 

and is known as the asymmetric a-model [||. For the symmetric case, a\ = a 2 , it is identical to the splitting function 
(^), except for the missing small anticorrelation; as a consequence the various multiplier distributions resulting from 
this symmetric case are more or less identical to those shown in Figs. 7 and 8. For the asymmetric case we choose two 
different sets of parameters: ct\ = 0.5, a 2 = 0.3 and a\ = 0.3, a 2 = 0.65. As depicted in Figs. 14a+b they represent 
two different asymmetries: the former enhances weights, which are larger than 1/2, whereas the latter enhances 
weights, which are smaller than 1/2. The two sets of parameters have been adjusted in such a way that the limiting 
form of the unconditional left multiplier distribution comes close to the experimentally observed parametrisation @ . 
This is illustrated in Figs. 15a+b. Also shown is the limiting form of the unconditional centred multiplier distribution; 
both parameter cases lead to a narrow distribution as central multiplier values around 1/2 are more enhanced than 
for left/right multipliers. Hence, as far as unconditional multiplier distributions are concerned the two different 
asymmetric versions for the splitting function lead to more or less indistinguishable results. For conditional multiplier 
distributions this is going to change! 

The limiting form of the conditional L/R multiplier distribution practically shows no dependence on the parent 
multiplier once the first set of parameters a\ = 0.5, a 2 = 0.3 is employed. This is exemplified in Fig. 15c. The second 
set of parameters «i = 0.3, a 2 = 0.65 reveals a different behaviour (see Fig. 15d): for a large i/i?-parent multiplier 
the conditional L/R multiplier distribution is significantly broader than for a small L/i?-parent multiplier; this trend, 
which would have also resulted without application of the scheme to restore translational invariance, is in agreement 
with experimental observations. - Also for the conditional centred multiplier distributions the two sets of parameters 
ot\, a 2 lead to different results: for a large/small parent multiplier the conditional centred multiplier distribution 
leads to a shift towards larger/smaller daughter multipliers. Whereas this positive correlation is observed for both 
parameter sets and is a consequence of the application of the scheme to restore translational invariance alone, only 
the second one leads to a further broadening/narrowing of the conditional centred multiplier distribution, as observed 
from data. 

These findings reveal that not every non-energy flux conserving splitting function qualifies to reproduce the various 
multiplier distributions. Moreover, in order to reproduce at least qualitatively the subtle effects observed in the 
conditional multiplier distributions, a certain asymmetry should be included in the splitting function underlying the 
discrete cascade process. This is in accordance with our intuitive picture of a three-dimensional cascade process 
observed in a one-dimensional world. In the next two Subsections we will test two popular models with asymmetric 
splitting functions: the log-normal ||l4] , |l5f and the log-Poisson model [p~6| p~8|| . 



B. log- normal weight distribution 

For the log-normal model the weight distribution entering the factorised splitting function is given by 



2e ff 1 
P(W L/R ) = -7=- exp <^ - —J 



XnW L/R + (\n2+^<7 2 ) 



(24) 



where the normalisation is such that (W®, R ) — 1 and (Wl/r) — 1/2. It is a little unphysical because, in principle, 
values Wl/_r. > 4 can be drawn; however, in practice the associated probability is negligible. The free parameter a is 
chosen such that the resulting limiting form of the unconditional left multiplier distribution, obtained after application 
of the scheme to restore translational invariance, comes close to the parametrisation (||); this gives a 2 = 0.45. Consult 
Fig. 16a, where also the resulting limiting form of the unconditional centred multiplier distribution is depicted; again 
the latter is more narrow than the former. 

Comparing the log-normal weight distribution (Fig. 14c) with the two asymmetric a-model weight distributions 
considered in the previous Subsection (Figs. 14a+b), the asymmetry of the log- normal distribution is more along 
the asymmetric a-model distribution with the second set of parameters a\ — 0.3, ot\ — 0.65 than the other set. If 
asymmetry towards smaller weights is indeed a criterion to reproduce the various experimentally observed conditional 
multiplier distributions on an at least qualitative level, then the log-normal model should also qualify. Figs. 16b+c 
illustrate the resulting limiting forms of the various conditional multiplier distributions. The conditioning on a 
left /right as well as centred parent multiplier leads to exactly the same qualitative effects as we have already realized 
for the "good" asymmetric parametrisation of the a-model. 
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C. log-Poisson weight distribution 

With the assumption that the most intermittent structures in three-dimensional turbulence are one-dimensional 
isolated vortex lines the scaling exponents r(n) = — Ci[n(l — C2) — (1 — c§)] with c\ = 2, C2 = 2/3 have been derived 
for the energy dissipation field ((e { ^) n ) = (L/2^) T ^ (T|. An inverse Laplace-transform directly leads to the weight 
distribution 

00 

P(W L/R ) = 2- Cl ^ Clln2rS ( Wl/r - 2Cl(1_C2Klc ") > ( 25 ) 

771=0 

which is of log-Poisson type 
function p(Wl,Wr) = p(Wl)p(Wr) for the evolution of a binary discrete cascade process. 

Observe, that there is no free parameter in the weight distribution ( p5| ) as long as parameters c\, C2 are fixed 
by general arguments. Hence, we can not finetune it to reproduce, for example, the unconditional left multiplier 
distribution. The more remarkable it is, that the respective limiting form obtained after application of the scheme 
to restore translational invariance comes close to the parametrisation (||); see Fig. 17a. Also the limiting form of the 
unconditional centred multiplier distribution is more narrow in comparison with the former distribution. 

Comparing the weight distribution ( p5| ) with the previously discussed weight distributions, see Fig. 14, it appears 
that the log-Poisson distribution clearly favours weights larger than 1/2 over smaller values. In this respect it more 
or less resembles the asymmetric binomial distribution with parameters a\ — 0.5, 0,2 = 0.3 and, hence, we expect 
no big effects in the conditional multiplier distributions. As illustrated in Figs. 17b+c this is indeed the case: the 
conditioning on a left/right parent multiplier leads only to relatively small deviations from the respective unconditional 
distribution. The conditional centred multiplier distributions again reveal a positive correlation between daughter 
and parent multiplier, but no additional broadening/narrowing is observed for a conditioning on large/small parent 
multiplier values. 

Instead of fixing the parameters c\, C2 by general arguments |jl6|| , they can also be used as free fit parameters. A 
comparison of scaling exponents for longitudinal velocity structure functions obtained from direct numerical simulation 
of low Reynolds number Navier-Stokes turbulence [|lj with those predicted from the log-Poisson cascade model 
suggests the parameter choice c\ = 9.3, C2 = (0.947) 3 = 0.849. With this choice the weight distribution ( p5| ) becomes 
more symmetric around Wl/r = 1/2 and, as a consequence, the various unconditional and conditional multiplier 
distributions nearly behave like those depicted in Figs. 7 and 8 for the symmetric binomial model with splitting 
function (JsJ) ; the multiplier correlations remain weaker than for the asymmetric binomial model (ot\ = 0.3, 012 = 0.65) 
or the log-normal model. 

VI. CONCLUSIONS 

Fully developed Navier-Stokes turbulence is a three-dimensional nonlinear process. With present experimental 
techniques it is not feasible to record this 3+1-dimensional spatio-temporal dynamics directly; standard experimental 
observations record time series of the velocity field in one point, which according to Taylors frozen flow hypothesis can 
be reinterpreted as a one-dimensional spatial cut through the three-dimensional velocity field at a given time. This 
reduction in dimensions is an important point to keep in mind when analysing and interpreting the "one-dimensional" 
data. - In this Paper we have discussed the implications of this dimensional reduction on one specific inertial-range 
observable, the so-called multiplier distributions associated to the energy dissipation field, by employing discrete mul- 
tiplicative branching processes. One-dimensional versions of the latter should not locally conserve the energy flux 
from large to small scales; moreover the splitting function should reflect an asymmetric weight distribution. As a 
consequence of this, the dynamical forward field used in the cascade evolution is not equal to the averaged backward 
energy dissipation field. This leads to deviations from perfect moment scaling [ p2[ . Unconditional multiplier distribu- 
tions obtained from the backward field do not match any longer the weight distributions employed during the forward 
evolution, but converge to scale-independent fix-point limiting forms which already come close to the experimentally 
deduced multiplier distributions. An even better agreement is obtained once also homogeneity is restored in the hi- 
erarchical cascade processes. Then, not only unconditional but also conditional multiplier distributions qualitatively 
match the experimentally observed distributions. We conclude that the correlations observed in various multiplier 
distributions are not in conflict with simple discrete multiplicative branching processes, but are a consequence of 
a one-dimensional homogenous observation of a three-dimensional non-homogeneous and strictly selfsimilar cascade 
process. 

12 



1 7 flq] and which we are now going to use as input into the bivariate factorised splitting 



So far the comparison between these simple cascade models and data had only been qualitatively. One natural 
question to ask now is: what is the correct, or at least best, one-dimensional splitting function in accordance with 
data? This study has already shown, that the asymmetric binomial model (with the second set of parameters) 
and the log-normal model yield almost identical output for the various multiplier distributions discussed so far. 
Hence, fine-tuning of different asymmetric parametrisations for the splitting function only makes sense once further 
observables are included; these could be for example multiplier or wavelet correlations [ p2|j23| and n-point correlations 
of the energy dissipation field itself |l3| |. Note, that in this last reference the inverse problem, how to reconstruct the 
splitting function from n-point correlations, has been solved in the pure discrete hierarchical cascade picture; however, 
the important aspect of restoring homogeneity destroys the beautiful formalism. - Only all these observables taken 
together are able to judge about the relevance and qualities of (discrete) cascade and other empirical models! 



ACKNOWLEDGMENTS 



One of us (B.J.) acknowledges support from the Alexander-von-Humboldt Stiftung. We thank Jiirgen Schmiegel 
for several fruitful discussions. 



13 



[1] U. Frisch, Turbulence (Cambridge University Press, Cambridge, 1995). 
[2] B. Mandelbrot, J. Fluid Mech. 62 (1974) 331. 

[3] U. Frisch, P.L. Sulem and M. Nelkin, J. Fluid Mech. 87 (1978) 719. 

[4] D. Schertzer and S. Lovejoy, in Turbulent Shear Flow 4, eds L.J.S. Bradbury et.al., (Springer, Berlin, 1985), p. 7. 

[5] C. Meneveau and K.R. Sreenivasan, Phys. Rev. Lett. 59 (1987) 1424. 

[6] C. Meneveau and K.R. Sreenivasan, J. Fluid Mech. 224 (1991) 429. 

[7] A.B. Chhabra and K.R. Sreenivasan, Phys. Rev. Lett. 68 (1992) 2762. 

[8] K.R. Sreenivasan and G. Stolovitzky, J. Stat. Phys. 78 (1995) 311. 

[9] G. Pedrizzetti, E.A. Novikov and A. A. Praskovsky, Phys. Rev. E 53 (1996) 475. 
[10] M. Nelkin and G. Stolovitzky, Phys. Rev. E 54 (1996) 5100. 
[11] M. Greiner, J. Giesemann, and P. Lipa, Phys. Rev. E 56 (1997) 4263. 
[12] E.A. Novikov, Appl. Math. Mech. 31 (1971) 231; Phys. Fluids A 2 (1990) 814. 

[13] M. Greiner, H. Eggers and P. Lipa, Phys. Rev. Lett. 80 (1998) 5333; M. Greiner, J. Schmiegel, F. Eickemeyer, P. Lipa and 

H. Eggers, Phys. Rev. E 58 (1998) 554. 
[14] A.M. Obukhov, J. Fluid Mech. 13 (1962) 77. 
[15] A.N. Kolmogorov, J. Fluid Mech. 13 (1962) 82. 
[16] Z.-S. She and E. Leveque, Phys. Rev. Lett. 72 (1994) 336. 
[17] B. Dubrulle, Phys. Rev. Lett. 73 (1994) 959. 
[18] Z.-S. She and E. Waymire, Phys. Rev. Lett. 74 (1995) 262. 
[19] B. Jouault, P. Lipa and M. Greiner, Phys. Rev. E, in press. 
[20] A. Bialas and R. Peschanski, Phys. Lett. B 207 (1988) 59. 
[21] S. Grossmann, D. Lohse and A. Reeh, Phys. Fluids 9 (1997) 3817. 
[22] M. Greiner, J. Giesemann, P. Lipa and P. Carruthers, Z. Phys. C 69 (1996) 305. 
[23] M. Greiner, P. Lipa and P. Carruthers, Phys. Rev. E 51 (1995) 1948. 



14 




FIG. 1. Schematic summary of the experimental findings with respect to the scale-invariant L/R- multiplier dis- 
tributions: (thick solid line) unconditional distribution p(Mj?^ R l ) = p(M^^\0 < M^ R < 1), (grey) distribution 
p{M^t R \Q < M_ < M^J R < M < 1/2) conditioned on a small L/i?-parent multiplier, and (hatched) distribution 
p(M^ / + 1) |l/2 < M < M^L < M < 1) conditioned on a large L/ii-parent multiplier. 
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FIG. 2. Schematic summary of the experimental findings with respect to the scale-invariant centred multiplier 
distributions: (thick solid line) unconditional distribution p(M^? +1 ') = p{M^ +r> |0 < Mq^ < 1), (grey) distribu- 
tion p(M^ +1) |0 < M < M^ < M < 1/2) conditioned on a small C-parent multiplier, and (hatched) distribution 
p(M^ +1) |l/2 < M < Mp' < M < 1) conditioned on a large C-parent multiplier. For comparison, the unconditional 
L/R- multiplier distribution p(M^^) of Fig. 1 is also shown as the thick dash-dotted curve. 
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FIG. 4. Cascade evolution of one energy dissipation field realisation from large to small scales (left column) vs. integrated 
energy dissipation field from small to large scales (right column). The non-energy flux conserving splitting function (H) has 
been employed with a = 0.4 and J = 10. Only every second scale j = 0, 2, . . . , 10 is shown. 
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Left multiplier distributions 
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Centred multiplier distributions 




M c (i) 



FIG. 5. Convergence of the (a) left and (b) centred unconditional multiplier distributions p(M^ ) and p(Mq ) of the 
one-dimensionally observed three-dimensional binomial model (^) to a quasi-continuous limiting form. Parameters are J = 9, 
a — 0.4. For comparison the experimentally deduced Beta-function parametrisation (^|) of the left (right) unconditional 
multiplier distribution with (3 = 3.2 is shown as a solid line. 
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Conditional multiplier distributions 




FIG. 6. Conditional (a) left/right and (b) centred multiplier distributions for scale j = 3 obtained from the splitting function 
(|) with a = 0.4 and J = 9. (a) p(M$ R \M < M^" 1 ' < M) for M = 0.2, M = 0.4 (grey) and M = 0.6, M = 0.8 (hatched); 

(b) p{M^\M < M^ _1) < M) for M = 0, M = 0.5 (grey) and M = 0.5, M = 1 (hatched). For comparison the experimentally 
deduced Beta-function parametrisation (^) of the left (right) unconditional multiplier distribution with (3 = 3.2 is shown as a 
solid curve in (a). 
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Left multiplier distributions 
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Centred multiplier distributions 




M c (i) 



FIG. 7. Same as Fig. 5, but now after restoration of translational invariance. 
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Unconditional multiplier distributions 
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FIG. 9. Convergence of the left (hatched) and central (grey) unconditional multiplier distributions p(M^ ) and p(M@ ) 
of the Beta-distribution cascade model ( |ll| ) to respective limiting forms. Parameters are J = 9, /3 = 3.2. For comparison the 
experimentally deduced Beta-function parametrisation (^) of the left (right) unconditional multiplier distribution with /3 = 3.2 
is shown as a solid line. 
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Conditional multiplier distributions 




FIG. 10. Conditional left (a) and centred (b) multiplier distributions for scale j = 3 obtained from the splitting function 
(0) with P = 3.2 and J = 9. (a) p(M^L|M < M^~^ < M) for M = 0.2, A? = 0.4 (grey) and M = 0.6, M = 0.8 (hatched); 

(b) p{M^\M < M^ _1) < M) for M = 0, M = 0.5 (grey) and M = 0.5, M = 1 (hatched). For comparison the experimentally 
deduced Beta-function parametrisation (^) of the left (right) unconditional multiplier distribution with (3 = 3.2 is shown as a 
solid curve in (a). 
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Unconditional multiplier distributions 




Conditional multiplier distributions 




FIG. 13. Conditional left (a,b) and centred (c) multiplier distributions for scale j = 3 obtained from the splitting function 
(0) with /3 = 3.2 and J = 9 after restoration of translational invariance. (a) p(M^j R \M < M { ^ } < M) for M = 0.2, M = 0.4 
(grey) and M = 0.6, M = 0.8 (hatched); (b) p(M^ R \M < M ( ^ R 1] < M) for M = 0, M = 0.2 (grey) and M = 0.8, M = 1 

(hatched); (c) p(M^\M < M^~ 1] < M) for M = 0, M = 0.5 (grey) and M_ = 0.5, M = 1 (hatched). For comparison the 
experimentally deduced Beta-function parametrisation (g) of the left (right) unconditional multiplier distribution with j3 = 3.2 
is shown as a solid curve in (a,b) . 
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FIG. 15. Various multiplier distributions resulting from the asymmetric binomial model after restoration of translational 
invariance. Left column is for model parameters Qi = 0.5, a?2 = 0.3, J = 9, and right column is for ai = 0.3, «2 = 0.65, J = 9. 
First line: unconditional left and centred multiplier distributions p(Ml) and p(Mc) at j — 3; second line: conditional left/right 
multiplier distribution p(M ( L j J R \M_ < M^" 1 ' < M) for M = 0.2, M = 0.4 (grey) and M_ = 0.6, M = 0.8 (hatched) at j = 3; 

third line: conditional centred multiplier distribution p{M^ ) \M_ < M^~ 1] < M) for M_ = 0, M = 0.5 (grey) and M = 0.5, 
M = 1 (hatched) at j = 3. For comparison the experimentally deduced Beta-function parametrisation (H) of the left (right) 
unconditional multiplier distribution with /3 = 3.2 is shown as a solid curve in each plot. 
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FIG. 16. Various multiplier distributions resulting from the log-normal-model after restoration of translational invariance; 
parameters are a 2 — 0.45, J = 9. (a) unconditional left and centred multiplier distributions p(Ml) and p(Mc) at j = 3; (b) 
conditional left/right multiplier distribution p{M^j R \M < M^~ R 1] < M) for M = 0.2, M = 0.4 (grey) and M = 0.6, M = 0.8 

(hatched) at j = 3; (c) conditional centred multiplier distribution p(M^ j) |M < M^ _1) < M) for M = 0, M = 0.5 (grey) and 
M = 0.5, M = 1 (hatched) at j = 3. For comparison the experimentally deduced Beta-function parametrisation (B) of the left 
(right) unconditional multiplier distribution with f3 = 3.2 is shown as a solid curve in each plot. 
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FIG. 17. Various multiplier distributions resulting from the log-Poisson-model after restoration of translational invariance; 
parameters are ci = 3c2 = 2, J = 9. (a) unconditional left and centred multiplier distributions p(Ml) and p(Mc) at j = 3; (b) 
conditional left/right multiplier distribution p(M^J R \M < M^~^ ] < M) for M = 0.2, M = 0.4 (grey) and M = 0.6, M = 0.8 

(hatched) at j = 3; (c) conditional centred multiplier distribution p{M^\M_ < M^ _1) < M) for M = 0, M = 0.5 (grey) and 
M = 0.5, M — 1 (hatched) at j = 3. For comparison the experimentally deduced Beta-function parametrisation (B) of the left 
(right) unconditional multiplier distribution with f3 = 3.2 is shown as a solid curve in each plot. 
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